Asymmetric long Josephson junction acting as a ratchet for a quantum field 
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We study the escape rate of flux quanta in a long Josephson junction having an asymmetric spatial 
inhomogeneous critical current density. We show that such a junction can behave as a quantum 
ratchet when it is driven by an ac current in the presence of a magnetic field. The rectification 
gives rise to an onset of the dc voltage Vdc across the junction. The usual approach of particle-like 
tunneling cannot describe this rectification, and a quantum field theory description is required. We 
also show that under definite conditions the rectification direction and, consequently Vdc, can change 
its sign when varying the temperature T near the crossover temperature T* between the quantum 
and classical regimes. 

PACS numbers: 05.40.-a, 74.50.+r, 85.25. Cp 
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Introduction. — Solid state devices with asymmetric 
periodic potentials (acting as ratchets) are attracting 
considerable interesiji. Ratchets can produce a direct 
current when a time-dependent force (deterministic or 
random) with zero mean is applied to it. This effect 
occurs for both: thermal hopping and quantum tunnel- 
ing. Moreover, for some special shapes of the ratchet 
potential, the current can change its sign for decreasing 
temperatures, when quantum tunneling becomes domi- 
nant over thermal-activated hoppingS.. There are many 
different realizations of ratchets, both in nature and in 
artificial nano-devices, such as cold atoms, colloidal mag- 
netic particles, single-molecule optomechanical devices, 
fluxons in superconductors, and many other systems (for 
reviews, see, e.g., Ref.[l|). All these systems are effectively 
described by a quantum or classical particle moving in a 
periodic asymmetric potential. 

Here we propose a completely new ratchet system, 
which is described by a quantum field. Namely, we 
consider a long Josephson junction (JJ) (the junction's 
length D is comparable with the Josephson penetration 
depth Aj) with a spatially modulated critical current 
density idx) driven by an adiabatically slow ac current 
J = Jq cos fit, where the amplitude Jq is smaller than the 
junction's critical current Jc- The gauge-invariant phase 
difference <fi{x) plays the role of the field variable. Al- 
though J never exceeds Jc, the flowing current can give 
rise to sudden changes in the phase difference, both due 
to the thermal hopping of the flux quanta through the 
potential barrier and quantum tunneling. The probabil- 
ity per unit time (escape rate) r( J) of these transitions 
attains maxima T± at J = ± Jq, and one can expect that 
r+ = r_. However, the application of an external dc 
magnetic field H to the junction having current inhomo- 
geneity asymmetric with respect to the x direction gives 
rise to the asymmetry in tunneling probability, r+ ^ r_ . 
This leads to the appearance of a dc voltage Vdc across 
the junction, since Vdc = (ft/2e)(^) cx {T+ — r_). 
We calculate the escape rate F using the well-known 



method of imaginary-time trajectories at finite tempera- 
ture^. According to this approach, F can be written as 
F = ^exp(— i3), where the prefactor A is of the order 
of the Josephson plasma frequency LUp, and B = Sp/h. 
Here, Sp is the action of the system calculated for a peri- 
odic imaginary-time trajectory, ip{T-\-h(3) = '^{t), with a 
period tq — h(3, (3 — l/T. The effect of quantum tunnel- 
ing in J Js and stacks of intrinsic J Js in high-Tc supercon- 
ductors is now studied intensively, experimentally^ and 
theoretically^, both due to its fundamental interest and 
the possibility to use these systems in future applications. 
A ratchet based on a SQUID consisting of two equal JJs 
in series, coupled in parallel to a third junction, was pro- 
posed in Ref. [y. However, the particle-like approach used 
in most theoretical considerations of quantum tunneling 
in JJs is not appropriate here, and a field-theoretical de- 
scription must be used. 

We stress that the system under consideration does 
not correspond to a ratchet in the usual sense since the 
junction's potential is the periodic sine-Gordon potential, 
U (tp) oc — cos if, which is symmetric in tp. The difference 
between F_|_ and F_. and, consequently, a voltage rectifi- 
cation occurs here due to both: the "parametric" depen- 
dence of the potential U on x, and the spatial dependence 
of the field (p{x). Thus, the coordinate dependence of 
the phase difference is crucial for the rectification, even 
for relatively short junctions, when D < Xj. Here we 
exploit the approach developed in Ref. 0, which we gen- 
eralize here to finite temperatures. This now allows us to 
calculate, in the same manner, the escape rate F both in 
the quantum and classical regimes. If the junction width 
D is large, then the exponent B of the escape rate F 
turns out also to be large. However it should be not too 
large, when F becomes an experimentally non-observable 
quantity. This last condition means that the current am- 
plitude Jo should be close to the critical current^. We 
focus on the exponent B because changes in B (not in 
the prefactor A) describes the main change in F when B 
is large. In addition, here we do not consider dissipation 
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in the system when calculating B. The effect of dissi- 
pation does not change qualitatively the results obtained 
here. 

Escape rate. — Wc consider a Josephson junction in the 
inline geometry shown in the inset of Fig. [TJ Two super- 
conducting bars overlap a length D in the .x-direction. 
The external magnetic field H is applied in the y- 
direction. Let us first consider the junction biased by 
a dc current J. The generalization to the adiabatically 
varying ac current J cos fit is evident, and the conditions 
of adiabaticity will be given below. The spatial inho- 
mogcneity of the critical current density ic{x) in the x- 
direction can be realized, e.g., by a spatial variation of the 
thickness of the insulating layer between superconduc- 
tors, or slightly changing the shape of the junction itself. 
We assume that Jc is homogeneous in the y-direction. 
The system under study can be described by an eff'ec- 
tive ID Lagrangian for the phase difference iy9(r, x). The 
quasi-classical imaginary-time action of this system is 
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(1) 



In Eqs. ([T]), the x-coordinatc is normalized by Aj, r by 
1/ujp, d = D/Xj, P = 1/T, Ej = fulX.,L/{2e) is the 
Josephson energy (L is the junction's length in the y- 
direction), the bar above ic means spatial averaging, and 

= + ^^0. (2) 



The phase (/^(t, x) satisfies the equation of motion 
d'^ip d'^p ... 



(3) 



with the periodicity condition lp{t + hujpP.x) ~ lp{t,x). 
The boundary conditions for the phase difference 
(/?(t, x) are defined by the jCs term in the action ([1]): 
dip/dx\x=±d/2 ~ ±//2 + /i, where / and h are the dimen- 
sionless current and external magnetic field, respectively: 
/ = J/{icXjL), and h = cH / {A-KicXj). 

In general, it is hard to find a periodic solution to the 
nonlinear Eq. ([3|), even numerically. However, as men- 
tioned above, we only need to find a solution when the 
current / w /c. In this case, we can use the approach de- 
veloped in Rcf . 0. Namely, we seek a solution of the form 
P>{t,x) ~ ipo{x) + 'ip{T,x), where ipoix) is the steady-state 
solution of Eq. ^ corresponding to an energy minimum. 
Since / is close to the critical current, the energy bar- 
rier between the neighboring energy minima, ipo{x) and 
(po{x) ± 27r, is small, |'0(a;,T)| <C ipoix), and we can ex- 
pand the action ([T]) in powers of ■0 up to ip^. Then, we 



represent '0(t, x) in the form of a scries 
3mo 
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where and ipnix) are the eigenvalues and orthogonal 
eigenfunctions of the operator— 



V = -d'^/dx^ + g{x) cos ipo{x) , Vipn = /^„Vn , (5) 



and 



I'd/2 

Unrnk = / dx g{x) SlTl p}o{x)'4'n{x)lprn{x)lpkix) . 
J -d/2 



(6) 



Expanding the equation of motion ([3]) in powers of "0, 
up to ip'^, with if) in the form (jl]), multiplying it by tpn{x), 
and integrating over x, we obtain a system of ordinary 
differential equations for the collective coordinates a„(77) 
of the field -;/', where we introduce the imaginary-time 
variable i] ~ ^/J^ot. It can be shown that due to the 
proximity of / to 1^ two conditions for the eigenvalues 
/i„ are possible: either <^ fin, n > for relatively 
short junctions {d < 1) or large fields {h > 1), or /ig ~ 
Ml ^ Mn, > 1 for long junctions {d> 1) a7T.d small fields 
{h <^ 1). Due to these inequalities, one can neglect all 
the equations in the system of equations for q;„, except 
the first two (for details, see Ref. 0)- Thus, the system 
of equations for and ai takes the form 



dV 
dao 



= 0. 
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where the dot means derivative over 77, and the potential 
V{aQ, ai) can be written as 
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Wc should find a periodic solution of the system ([7]) with 
period r/o = y/jM^fujp/T ("bounce" solutions). The expo- 
nent B = Sp/h of the escape rate F — Acxp{—B) can be 
expressed through the functions Oii{rf), z = 0, 1 as: 



000 j 



5]^-y(ao,ai) 



,A=^. (9) 



Thus, we reduce the field-theory problem to the prob- 
lem of one particle moving in the potential V{aQ,ai) 
in two dimensions, where the a^'s play the role of the 
particle's generalized coordinates. Let us first consider 
the case of T ^ 0, when quantum tunneling prevails. 
Our analysis shows that when d > dc{I,h) ^ 4, there 



exist three bounce solutions to the system (I7|): 



s = 0,±1. In analogy to the case 7(x) = /i = con- 
sidered in Ref. 0, one can say that the solution a[p^^\ri) 

(q!j- ^\r])) corresponds to the formation of a vortex (an- 
tivortcx) nucleus at the left (right) junction edge, while 
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FIG. 1: (Color online) The magnetic field dependence of the 
escape rate's exponent B for fixed current ratios I/Ic(h,d) 
in the uniform junction [7(2;) = 0], calculated according to 
Eq. @; d = 1.8 and T = 0. The inset shows the schematic 
geometry of the Josephson. 

the solution a^"' (77) describes the tunneling of as a 
whole (for details, see Ref. 0). When ^{x) = h = 0, 
the solutions have the following symmetric properties: 

= Oj '^i~^\v) = ~'^i^^\v)- Thus we have 
three channels of tunneling, with probabiUties F^*' cx 
exp(-B('')), with B^-i) = B'-+^\ and the total proba- 
bility becomes F = F''*^ The applied magnetic field 
breaks the vortex-antivortex symmetry [B^^^^ = S'+^^J 
making one of these channels more favorable. However, 
the total escape rate F is still symmetric with respect 
to the direction of the dc current if 7(x) ~ 0. For spa- 
tially inhomogeneous junctions with 7(x) 7^ 0, we have 
F(-|-J) ^ r(— /), and the rectification arises. Since the 
main contribution to the total F comes from the term 
corresponding to the minimum of B^'^\ we will assume 
below that F cx cxp(— i?), where B ~ min (B(^'). When 
d < dc{I, h) there is only one solution to the system ([7]), 
a|°^(?7). When h^Q and/or 7(0;) ^ 0, we have 7^ 0, 
and. in contrast to the case studied in Ref. 0, a 2D con- 
sideration is required here. 

All these properties of the bounce solutions af^ survive 
at finite temperatures up to some value T*, which is the 
crossover temperature between the quantum and classical 
regimes. At low temperatures, T < T* , the exponent 
B only slightly decreases with T (thermally-stimulated 
tunneling). When T > T*, there arc only imaginary- 
time independent solutions to Eq. ([7]), a['^\r]) = a['^\ 

where the points {alj^\ a[^^) correspond to the extremes 
(minimum or saddle-points) of the potential V{ao,ai) 
(here s takes the values s = 0, ±1, if d > dc{I,h), and 
s — otherwise). Thus, when T > T* the exponent B 
can be written as 



B = 



9A/ii] hujp\Vr\ 



T 



Vr 



(10) 



The crossover temperature T* is the temperature where 



FIG. 2: (Color online) The magnetic field dependence of 
B{+I, h) (red solid curve) and B[—I, h) (blue dashed curve) 
for a fixed ratio I /Ic ~ 0.96, calculated for 7(1) = 0.014a;. 
The relative ratio 1 — B{—I, h)/B{+I, h) as a function of h is 
shown by the green d otted- dashed curve. Here the parameters 
are d = 2.5, T = 0, {^f(x)y^^= 0.01. 



the period of the bounce solution, 770 = ^/jM)fTi^p/T , be- 
comes equal to the period of infinitesimal oscillations 
near the extreme point (ao'^\ ol^i^)^ corresponding to the 
maximum of V {6l'^\ ol"-^^ ) . The latter one is equal to 
27r/V'A^, where is the positive^^ eigenvalue of the 

matrix d'^V / daidaj calculated at the point (q:q'^\ Q^i'^^)- 
As a result, the crossover temperature can be written as 
T* ~ hujp\/ ^o\t /{2n). 

All results above were obtained for a dc current. We 
can also use all the above formulas for a slowly-varying 
current / = /q cos Jli, assuming that at any time t the 
current can be considered as a dc one. This adiabatic 
limit is valid when the frequency Q is much smaller than 
the inverse imaginary-time period, I/tq = T/h. There- 
fore, we obtain the condition of adiabaticityii: hil <^ T. 

Results and discussions. — Let us first neglect the spa- 
tial distribution of the critical current density, assuming 
that 7(2;) = 0. In this case, the following condition is 
met B{—I, h) — B{I, h), and there is no rectification. In 
Fig. [U we plot the magnetic field dependence of the ex- 
ponent B {T = 0) for three fixed ratios I/Ic{h, d), where 
the critical current Idh, d) depends on the magnetic field. 
It goes to zero when h achieves the critical field hc{d). 
Note the non-monotonic dependence of B on h, which 
is related to the change in characteristic properties of 
the static solution ipa{x) with growing h. The non-zero 
7(2:) breaks down the symmetry of B with respect to 
the direction of the current. In this case, we have two 
critical currents I^{h,d) and two critical fields h^{d), 
corresponding to positive and negative currents. Assum- 
ing that |7(x)| ^ 1, we calculate the exponent B in first 
order perturbation theory with respect to 7. Figure [2] 
shows the dependence of B{I, h) and B{—I, h) on /i, cal- 
culated for linear 7(x) at the fixed ratio I/Ic{h, d), where 
Ic{h,d) is the critical current calculated to zeroth-order 
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FIG. 3: (Color online) The escape rate exponents B{+I) (red 
solid curves) and B{—I) (blue dashed curves) versus temper- 
ature for two different spatial 7(2:) distributions (shown in the 
inset): -yix) = 0.014a: (a) and 7(0;) = 0.17a; - 0.27a:^ (b). The 
differences B{+I) — B{~I) are shown by green dotted-dashed 
curves. The other parameters, d = 2.5, h = 2.25, I/Ic = 0.95, 
are the same for both panels. 

in 7. The degree of rectification, 1 — [B{—I, h)/B{I, h)], 
shown in Fig. [2l monotonically increases up to /i = , 
where I~ = 0. 

Let us now consider the temperature dependence of the 
escape rate exponent B. The curves B± — B{ztl) ver- 
sus T, calculated for two different spatial profiles 7(x), 
are shown in Fig. [31 The difference — which 
defines a ratchet effect as a function of temperature is 



also shown in Fig. [3] For linear 7(1) (Fig. [3^) this dif- 
ference is always positive, while for the more asymmetric 
profile 7(2;) = 0.17x — 0.27x^ (which is shown in the in- 
set to Fig. [3]d) it changes sign at temperatures near the 
crossover temperature between the quantum and classi- 
cal regimes. There is a big difference in absolute values 
of the ratchet effect between these two cases: the relative 
ratio (B+/B- - 1) for 7(0;) = 0.17x - 0.27x^ (when it 
changes sign) is much smaller than for linear ^{x). It can 
be enhanced by optimizing the profile for 7(0;). The effect 
of the profile idx) on the rectification direction can be 
understood following the explanation given in Ref. [2: for 
thermoactivated hopping only the height of the poten- 
tial barrier is essential, while the probability of quantum 
tunneling depends also on the thickness of this barrier. 

The sign of the difference {B+ — B-) defines the sign of 
the rectified voltage Vdc- Thus, the effect of the change of 
sign of (i3_|_ — i?_ ) can be observed by measuring the dc 
voltage as a function of temperature. The value (5+ — 
B^) can be calculated using data of the escape rate F 
obtained cither in ac or dc current measurements. 

Conclusions. — We have proposed a new type of 
ratchet system based on a long Joscphson junction with 
a spatially inhomogeneous critical current density ic{x). 
The exponent of the escape rate of the phase difference ip 
was calculated as a function of: temperature, dc magnetic 
field, and dc or slow ac current. We have shown that due 
to both the magnetic field and the spatial inhomogencity 
of ic, the escape rate becomes asymmetric with respect to 
the direction of the current. This leads, in particular, to 
the appearance of a dc voltage when the system is biased 
by an ac current. We have also shown that, for definite 
shapes of ic{x), the rectified voltage changes sign at a 
temperature near the crossover temperature between the 
quantum and classical regimes. 
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® In real time, ipn correspond to linear standing-waves with 
frequencies LOn ~ ^/JI^Up excited over the solution ifo{x). 
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and d < dc, so the matrix of second derivatives of V at this 
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When T < T* , one can use a weaker inequality, hfl <^ T* , 
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